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1. INTRODUCTION 

Consider an animal that is facing a binary decision (YES/NO, 
GO/NO-GO) in response to the presence/absence of a stimulus 
(light/no light, sound/no sound). The animal makes its decision 
based on the firing pattern of its sensory neurons. What is the 
best decision strategy, one that minimizes the response time whUe 
keeping the number of decision errors within an acceptable pre- 
defined limit? In a wide set of situations the best decision strategy 
is well known and it is based on Wald's sequential probability ratio 
test (SPRT) (Wald, 1945; Wald and Wolfowitz, 1948). Models of 
neuron-based decision that are based on SPRT have been suc- 
cessful at predicting the response time histograms of animals 
that are engaged in decision tasks, and neurophysiological cor- 
relates of key stages of the computation have been found (Gold 
and Shadlen, 2007), suggesting that the brain may be using this 
decision strategy. Here I explore the predictions of this theory 
when very few, one in the limit, neurons are engaged in the deci- 
sion process. I find that response time histograms become highly 
structured, something that ought to be easy to reveal with simple 
behavioral/psychophysical experiments and might provide addi- 
tional insight into how the brain computes decisions from sensory 
inputs. 

In order to make this paper self-contained, in Section 2, I 
review SPRT and in Section 3, 1 review how to compute such deci- 
sion based on the firing pattern of Poisson (or Bernoulli) neurons. 
In Section 4, I explore an example with physiologically plausible 
constants. In Section 5, 1 compute the properties of response time 
histograms for a single neuron. In Section I briefly discuss the case 
of multiple neurons involved in a decision. Section 8 describes 
the computational experiments. I conclude in Section 9 with a 
discussion of the main observations. 

2. OPTIMAL BAYESIAN OBSERVER 

I assume that the animal is sensing, and responding to, an under- 
lying binary state of the world, e.g., the presence/absence of a 



The histograms of response times of optimal YES/NO decisions that are computed from 
a single sensory Poisson neuron are highly structured. In particular, response times in NO 
decisions are quantized to a small set of times, while response times in YES decisions 
have a multimodal structure. Both the times of NO decisions, as well as the modes of the 
the histogram of YES decisions, are associated to the number of action potentials that 
were necessary to reach the decision. Their value is a function of the firing rate of the 
neuron in response to the states of the stimulus. 

Keywords: decision, ideal observer model, poisson neurons, response time distributions, response time modeling, 
sequential probability ratio test, SPRT 



given stimulus. I will call C = 0 one state (stimulus absence) 
and C = 1 the other state (stimulus presence). The animal's 
YES/NO decision is computed from the stimulus. The animal is 
rewarded if the response is appropriate, i.e., (C = 0) — > NO and 
(C = 1) — ^-YES, and the animal is not rewarded otherwise. 

Call Xt the measurements on which the decision is based. I will 
assume that Xt is the firing pattern of one sensory neuron which 
will be firing vigorously if the stimulus is present, and will fire at a 
low level of activity if the target is absent. Here f is a discrete index 
of time, with observations starting at f = 1, thus, if Xt is the piece 
of information that is acquired at time t, then Xt = {xi, . . . , Xt}. 
As I will explain below, the theory may be developed both in 
discrete time and in continuous time; the two formalisms yield 
identical results. 

Wald's SPRT strategy works by repeating at each time t the 
following two steps: 

First, compute the log ratio of the probability that C = 1 is 
true, vs the probability that C = 0 is true, given the available 
dataX,: 

P(C=l\X,) 

Rt = log (1) 

^P(C = 0|X() 

Second, compare Rt to two thresholds, tq and ti (typically tq < 
0 < Ti). If Rt > Ti, then C = 1 is much more likely than C = 0, 
and decision YES is made. If instead Rt < tq, then decision NO is 
made. The third possibility is that tq < Rt < ti, i.e., the ratio is 
between thresholds. In this case no decision may be made because 
the information is insufficient, and one waits for the next piece of 
evidence Xt+i. The process is repeated until a decision is made. 
In summary: 



repeat for each t = I. 
obtain X(, compute Rt, if 



until decision reached: 

Rt > Ti — > decide YES 
Rt < To — > decide NO 
Rt e (To, Ti) — > wait for Xt+i 
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FIGURE 1 1 Probability of error vs. threshold on the log-likelihood ratio 
Rt (see Equations 2, 3). 
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FIGURE 2 I Bernoulli representation of a sequence of action potentials. 

In tlie black plot (above) the action potentials produced by a neuron over 
time are indicated with vertical bars. The binary X; € {0, 1 ) variable (blue 
plot, below) is obtained by dividing time into (short) bins of length At; 
Xf = 1 if an action potential is observed in the corresponding time interval f, 
and it is equal to 0 otherwise. The bins should be short enough that the 
probability of observing multiple action potentials in one bin is very small 
(e.g., At = 1 ms). By taking the limit At — > 0 one obtains a Poisson model 
of neuronal firing. 



Wald proved that this strategy is optimal (Wald, 1945; Wald and 
Wolfowitz, 1948), i.e., no strategy with the same error rate yields 
faster decisions. A review of this strategy in the context of neu- 
rophysiological models of decision may be found in Gold and 
Shadlen (2007). 

The two thresholds tq , ti determine the error rates, and may 
be computed from the expression for Rt in Equation 1: 



Rt = log 



P(C=1\X,) 
P(C = 0\Xt) 



log 



P(C=l\Xt) 
l-P(C=l\Xt) 



Table 1 | Parameters used in the example of section 4 and in the 
simulations shown in Figures 4, 5. For convenience, base 10 was 
used for both exp and log. 

Independent variables: 

= 1 spike/s 
= 10 spikes/s 



P(C = 0) = 0.5 
P(C= 1) = 0.5 



Firing frequency of neuron 
when C = 0 

Firing frequency of neuron 
when C = 1 

Prior probability of C = 0 
Prior probability of C = 1 



Derived variables: 

ro = log(P(C=1)/P(C = 0)) = 0 
r(1 ) = log(k^ /Xq) = 1 

« = logio(e)(;Vo-^i)«^-3.9s-i 



Initial value Ro of the 
diffusion (Equation 4) 
Diffusion increment, action 
potential at time t 
(Equation 7) 

Diffusion drift when no action 
potential is observed 
(Equation 6) 




(2) 



FIGURE 3 I Schematics of a diffusion resulting from SPRT. The vertical 
axis indicates the log probability ratio R(. The diffusion is indicated in blue. 
In the absence of action potentials the diffusion drifts by & units/s and 
reaches the lower threshold at time RT^ (see section 5 and Equation 8), 
where a NO decision is made. When an action potential occurs the 
diffusion jumps by r(^). When one of the two thresholds t is crossed, a 
decision is made. In this case the decision is NO; it is indicated by a red 
dot. It is made at time RT!^ (see Equation 8). Notice that if the time of the 
action potential changes within the interval t e (0, RT^), the decision is 
made at the same time. A number of diffusions produced by a simulation 
with the parameters shown in Table 1 are shown in Figures 4, 5. 



where the denominator of the second fraction depends on the fact 
that C may only take two values, 0 and 1. Solving this expression 
forP(C= l|Xt) yields: 

P{C = 1 \Xt) = — / = g(R,) (3) 

1 + exp(-J?t(X()) 



Where g(-) denotes the logistic function. Thus, say, if = ti = 
2, that means that P{C = l\X) ^ 0.99 (Equation 3 and Figure 1), 
i.e., the false accept error rate is 1%. Conversely, if a 1% false 
accept error rate is deemed acceptable, then (Equation 2) maybe 
used to compute the upper threshold on i^j: ti ~ 2 . Similarly, 
if a false reject rate of 1% is reasonable, then tq = —2 (for 
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FIGURE 4 I Patterns of action potentials and ensuing decisions when NO decisions by red dots. One may notice one false alarm error in the 

C = 0. Action potentials are represented by black vertical bars along the last row. The parameters of the simulation are shown in Table 1. For 

time axis. The corresponding diffusion RtiX) (Equation 4) is represented convenience, examples are sorted by the number of action potentials 

by a light blue line. When Rt crosses one of the two thresholds r a that were sufficient to reach a decision. For a frequency of each such 

decision is made. YES decisions are represented by green dots and event see Figure 6. 



convenience of mental calculation I am using base 10 for both where 
exp and log). Figure 1 shows the relationship of error probability 
and threshold magnitude. 

Computing Ri, Ri, Rt is easy from the class conditional 
probabilities P(Xt\C = 1), P{Xt\C = 0) if the measurements Xt 
are independent when the state of the world C is known. 
Using Bayes' theorem: 



Rt = log 



P(C= l\xu 



P(C = 0\xi,...,Xt) 

P(Xs\C=l) 

log 

5=1 



Xt) P(C = 1) 

— = loe 

^ P(C = 0) 



P(Xs\C=0) 
'0 + ^ r(Xs) = Rt^i+ r{xt) 



(4) 



:l0! 



r(xt) = log 



, PjC = 1) 
' P(C = 0) 
P(xt\C=l) 
P(xt\C=0) 
iiven evidence Xt 



(log ratio of priors) 

(log ratio of probabilities 



s= 1 



Equation 4 shows that the computation is a diffusion, i.e., it may 
be carried out recursively by updating the previous value Rt-i 
with a term r(Xf) that depends only on the current observation 
Xt, rather than on the whole set of observations Xt. The fact 
that some binary decisions may be implemented by the brain 
with a diffusion was first suggested as a phenomenological model 
(Ratcliff and Hacker, 1981) and later shown to be optimal, under 
appropriate conditions (Gold and Shadlen, 2001). 



Frontiers in Computational Neuroscience 



www.frontlersin.org 



April 2014 I Volume 8 | Article 42 | 3 



Perona 



Quantized response times 



RT=384 ms - 0 spikes 



RT=384 ms - 0 spil<es 



RT=640 ms - 1 spikes 



RT=640 ms - 1 spikes 






2 






ift 1 1 III III 


1 


K II 1 II nil 




r 


0 








-2 





Ml III 



INI 



0 0.5 1 1.5 2 

RT=27 ms - 2 spikes 



0 0.5 1 1.5 2 

RT=28 ms - 2 spikes 



0 0.5 1 1.5 2 

RT=199 ms - 3 spikes 



0 0.5 1 1.5 2 

RT=279 ms - 3 spikes 




III III 



Mil II 



II MM 



0 0.5 1 1.5 2 0 0.5 1 1.5 2 0 0.5 1 1.5 2 0 0.5 1 1.5 2 



RT=568 ms - 4 spikes 



RT=412 ms - 4 spikes 



RT=765 ms - 5 spikes 



RT=874 ms - 5 spikes 



2 
1 

-1 

-2 





2 






r 

|\jMi III 1 II 1 II 


1 
0 
-1 


kl 


1 II 1 1 lllll 1 




\|' 




-2 







2 






1 


• 


II II 1 II III II II 








0 


^^1 III 1 III 1 III 




-2 





0 0.5 1 1.5 2 0 0.5 1 1.5 2 0 0.5 1 1.5 2 0 0.5 1 1.5 2 

Time (s) 

FIGURE 5 I Patterns of action potentials and ensuing decisions when C = 1. Same parameters as in Figure 4. One may notice false reject errors in tlie 
first row. 



3. DECISIONS BASED ON ONE POISSON NEURON 

As soon as one knows the statistics P(Xt\C) of neuronal responses 
Xt to stimuli C one is able to compute the log likelihood ratio Rt 
explicitly. Here I will derive Rt assuming that neurons produce 
patterns of action potentials that are distributed with Poisson 
statistics, a model that has been shown to be useful in many 
instances (Seung and Sompolinsky, 1993; Jazayeri and Movshon, 
2006; Graf et al., 2011). More general Poisson-like models (Beck 
et al, 2008) may be used in the following analysis. They have the 
disadvantage of being more complex to implement and analyze 
and therefore I will confine myself to Poisson here. 

The equations are particularly simple if one uses a Bernoulli 
approximation to the Poisson distribution: call Xt the firing pat- 
tern of a neuron. I will assume that time has been discretized 
in small non-overlapping identical bins, e.g., time bins that are 
At = 1 ms long, and that the random variable Xt G {0, 1} repre- 
sents whether an action potential is observed during the interval 
corresponding to time -bin f, in which case Xt = I, otherwise 
Xf = 0 (see Figure 2). 



Consider a neuron whose expected firing rate is action 
potentials per second when C = 0, and action potentials 
per second when C = 1 (if one models the firing rate with 
a Poisson process, k is the expectation). In this case one 
may compute the value of r(xt) to be used in the diffusion 
Equation 4: 



r(xt) = log 



P(xt\C=l) 



•P(,xt\c = 0) 
the variable Xt takes only two values, 0 and 1. Thus: 

I -Pi 



r(0) = log 
r(l) = log 



1 -po 

Pi_ 
Po 



(5) 



where po and pi are the probabilities of detecting an action poten- 
tial in a given time bin when the world is in state C = 0 and C = 1 
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FIGURE 6 I Histograms of 2 ■ 10^ decision times for diffusions with 
parameters as in Table 1 and Figures 4, 5. The symbol C indicates tine 
stimulus and the symbol C indicates the decision; thus, the combination 
C = 0, C = 1 denotes false alarms and C = 1, C = 0 indicates false 
rejects. (Top-Left, Top-Center) The histogram of response times is sparse 
for C = 0, in line with the theoretical prediction (Section 5). In this case, 
decisions are more likely to be made with 0 action potentials, slightly less 



likely with 1, then 2 etc. (Bottom-left, Bottom-Center) The histograms for 
C = 1 is not sparse; however, it is lumpy. (Top-Right, Bottom-Right) The 

structure of the C = 0 and C = 1 histograms is clear when separate 
histograms are drawn for decisions made after k=0,...,6 action potentials. 
Notice that for C = 1 no such decisions involve fewer than 2 action 
potential, since 1 action potential is insufficient to cross the ti = 1.5 
threshold as discussed in Section 4. 



respectively. Thus, indicating with At the duration of a time bin, 
chosen so that kAt <^ 1: 

po ^ XoAt 
pi ~ XiAt 

The limit for Af — > 0 yields the continuous model. There is only 
one delicate point: in the limit r(0) — > 0; however, S = r(0)/ Af, 
the linear drift rate of the diffusion when no action potentials are 
observed, is different from zero. Thus: 

S = lim = lim — log — 

At^o At At^o Af I - po 

1 1-XiAf 

= = lo8(e)(^o - ^i) (6) 

Af^o Af 1 - ,>.oAf 



r(l) = log — = log — (independent of Af, thus the limit 

is trivial) (7) 

By taking the limit for Af — > 0 we obtained the exact equations 
for the Poisson model (Jazayeri and Movshon, 2006; Chen et al., 
2011). In Equation 6 the limit may be computed by consider- 
ing the Taylor expansion log(l +x)/log(e) = x + O(x^), where 
0(x^) indicates higher order terms that vanish when one takes 
the limit. When using base 10 for the log, as in Table 1 and in 
Section 4, then logiQ(e) ~ 0.4343. When using natural logs, then, 
of course, log(e) = 1. 

According to Equations 6, 7 the diffusion will mostly drift lin- 
early at a rate S, and present jumps of height r(l) whenever an 
action potential is observed (see Figure 3). 
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FIGURE 7 I When ki is reduced the spacing of the peal<s in the response time histograms becomes larger. Decisions take longer, since tine neuron fires 
less vigorously, and firing rates between C = 0 and C = 1 are more similar. Here Xo = =3; compare with Figure 6 where Xo = =10. 



Notice that if one assumes, without loss of generality, that X i > 
Xo, then the diffusion jumps will always be positive, i.e., upwards 
and the drift will always be negative, i.e., downwards, i.e., action 
potentials always contribute evidence toward C = 1, while quiet 
periods always contribute evidence toward C = 0. 

4. A CONCRETE EXAMPLE 

In order to develop one's intuition it is useful to explore a 
concrete case. Suppose that a neuron responds to stimuli with 
parameters indicated in Table 1, which correspond to a neu- 
ron firing briskly at 10 Hz when the stimulus is present (C = 
1) and at some "resting" level of 1 Hz when the stimulus is 
absent (C = 0). The diffusion starts at a value of Rq = 0, it is 
incremented by r(l) = 1 units whenever an action potential is 
observed (this is a rare event) and it drifts by 5 = — 9logjQ(e) ~ 
—3.9 units/s during the intervals where no action potentials 
are observed (see Figures). Simulations of this are shown in 
Figures 4, 5. 

Simple calculations allow one to estimate the expected 
response time for decisions when C = 1 and when C = 0. 



Suppose that we consider an error rate in the neighborhood 
of 3% acceptable, then we should use thresholds ti = 1.5 and 
To = —1.5 (Equation 3, Figure 1). Since one action potential 
increases the diffusion value by one unit, and 1 < ri, it will 
take at least two action potentials to reach a YES decision, while 
a sufficiently long time interval with no action potentials will 
lead to a NO decision (see Figure 3). We explore in the follow- 
ing the response time for YES and NO decisions when C = 1 
and when C = 0. A more systematic analysis is presented in 
Section 5. 

YES decision after 2 action potentials (C = 1) — Let's con- 
sider first the case where the decision is made after two action 
potentials. After two action potentials the value of the diffusion 
is Rt = 2 + St (f is the time that has passed since the begin- 
ning of observation). Thus a YES decision may be taken only 
if f2 < (2 — ri)/|3| ~ 128ms. The neuron fires with frequency 
Xi = W spikes/s, which makes it somewhat unlikely that two 
action potentials will be observed in 128 ms (the probability is 
about 0.23, which may be computed considering the Poisson dis- 
tribution with expectation X = fiAi). Thus, YES decisions that 
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FIGURE 8 I Diffusions for a decision involving two neurons. The neurons 
are symmetric witii respect to C: one fires more vigorousiy for C = 0 
(A.0 = 10, A.1 = 1 ) and tine otiier for C = 1 (io = 1 , = 10). In this case there 
is no drift (the two drift terms cancei) and the diffusion oniy moves when an 
action potential is observed. (Top row) examples of diffusions where a 
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decision is reached quickly, (IVIiddle row) examples of diffusions where 
decisions are reached after multiple action potentials. YES decisions are 
represented by green dots and NO decisions by red dots. (Bottom row) 
Histograms of response times. Notice that the response time histograms are 
neither sparse nor lumpy (compare with Figure 6). 



are based on two action potentials will only happen in a minority 
of cases. 

YES decision after 3 action potentials (C = 1) — Consider 
now the case where the YES decision is made after three 
action potentials. This decision may not be made when t < 
128 ms since in that case it would be made after the first 
two action potentials (see previous paragraph). Thus the deci- 
sion time is at least 128 ms. Furthermore, following the rea- 
soning in the previous paragraph, the decision time will at 
most be f3 = (3 — ri)/|S| ~ 384 ms. The probability that at least 
three action potentials are observed in less than 384 ms is 
about 0.47. 



NO decisions (C = 0) — Now consider the C = 0 case. If 
no action potentials are observed it will take tq/S ~ 1.5/3.9 ~ 
0.384 ms for a negative decision to be made. However, if one 
action potential is observed during this time, then the diffusion 
will be incremented by 1 and one will have to wait ~ 1/3.9 ~ 
0.256 s longer for a negative decision. When C = 0 it is quite 
likely that no action potentials are observed during a 0.384 s time 
interval since Xq = 1 s/s; the probability of this event is 0.68 
(computed from a Poisson distribution P(k = 0; X) with X = 
0.384). Thus, one would expect the majority of negative decisions 
to be made after no action potentials are observed and a time of 
about 0.384 s has elapsed. 
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FIGURE 9 I Decision times are shorter whien X^ax >> ^mm ■ (Left) Histogram of the response times for the two-neuron system explored in Figure 8. (Right) 

Median response time as a function of when = 1 spike per second. 



I simulated such a process and obtained histograms of decision 
times (Figure 6). The histogram of decision times for C = 0 is 
sparse. Only a few discrete decision times are observed. This issue 
is explored in Section 5. 

5. SPARSE RT HISTOGRAMS 

Each decision is made once a small number of action poten- 
tials has been observed. Without loss of generality, let's assume 
that the neuron fires more vigorously when C = 1 than 
when C = 0, i.e., ki > Xq =^ pi > pQ. Then r(l) > 0 > 5 (see 
Equations 6, 7). 

NO decisions — Since 5 < 0, decisions made after observing 
k = Q action potentials must be of type NO and are made in 
exactly RT^ = tq/S s (Figure 3). For = 1 a NO decision takes 
longer because the diffusion is incremented by r(l) when the 
action potential is observed (Figure 3); therefore, it takes — r(l) /S 
longer to reach the threshold, thus RT^ = RT^ - r(l)/S. The 
general expression for the NO decision time after k action poten- 
tials is therefore: 



Rn 



To 

J 



r(l) 



(8) 



Therefore, one would predict that when only one neuron is 
involved, decision times are sparse and no other decision time 
may be observed. 

Using the constants of Table 1 as in the example above yields: 
RTj^ = (1.5 - A:)/(9logio(e))), i.e., one would predict the follow- 
ing discrete decision times 384, 640, 895, 1151 ms etc. for k = 
0, 1, 2, 3, . . .. This is precisely what is observed in the simulation 
shown in Figure 6 (left). 

YES decisions — As shown in Figure 6, the histogram of 
response times is lumpy. It is possible to predict this observa- 
tion, and to see that each mode of the histogram corresponds 
to a different number k of action potentials. The time it takes 
to compute a YES decision that is based on k action potentials 
has a lower and an upper bound. Let's call the bounds RT^. ' and 



Y u 

RTj^ ' . The maximum time that it may take for a YES decision 
to be made is easy to compute: k action potentials increase the 
diffusion by fcr(l) units. To obtain a YES decision after k action 

Y u 

potentials in the amount of time RTj^ ' the diffusion will reach a 

Y u 

value equal to ti after RT^. ' seconds, which implies a downwards 

drift of aj^rf 



yields: 



fcr(l) — ri units. Solving for the upper bound 



rt: 



Ti-fc-r(l) 



(9) 



The lower bound may also be computed considering the fact 
that it is achieved when a decision based on — 1 action poten- 
tials is missed by a hair's breath and is followed immediately 
by another action potential, which overshoots the threshold. 
Thus: 



RT, 



Y, ri - (fc - 1) • r(l) 



(10) 



Using the constants of Table 1 yields the following predic- 



tions: RT2' ' 



: 0 s, RT, 



0.128 s, RT 



Y. I 



0.128 s, rt: 



0.384 s, J^Tj^- ' = 0.384, RT^^ 
shown in Figure 6 (right) are consistent with these predictions. 

If one compares equations 8 with 9 one can readily see that 
when Ti = — To the sparse RT for NO decisions and the "zeros" 
of the RT histogram of YES decisions are the same, i.e., YES 
decisions mostly happen at times that are different from NO 
decisions. 

One last question one may ask is how the various parame- 
ters (t and X) affect the spacing between the response times, 
and the height of the peaks, in the NO decision sparse his- 
togram. It is intuitive that the larger the spacing and the higher 
the second peak w.r. to the first one, the easier it will be to 
observe the sparse nature of the response time histogram in an 
in-vivo experiment where multiple sources of noise and vari- 
ability will tend to blur away the sparse/lumpy nature of the 



rY, U 



0.640 s etc. The simulation 
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histograms. From Equation 8 the spacing between peaks is 
equal to: 

^ _ r(l) _ log(Ai) - log(Ao) 
8 log(e)(Ai -Ao) 

therefore the spacing is maximum for ki — > Xq. An experi- 
ment where C = 0 does not stimulate the neuron, and C = 
1 stimulates the neuron lightly may therefore be the best 
option for revealing sparse response time histograms (see 
Figure 7). 

6. DECISIONS INVOLVING MULTIPLE NEURONS 

Decisions involving multiple neurons firing independently may 
be made using SPRT and Equation 4 (detailed equations for 
Poisson neurons may be found in Chen et al., 2011). Each 
action potential contributes to the diffusion independently of 
the other action potentials, and the contribution of each action 
potential is r'(l) = log(Aj) — log(/lg) (same as the case for 
a single neuron, see Equation 7), where i is the index of 
the neuron that generated the action potential. Thus, neu- 
rons that respond very differently to C = 0 and C = 1 (and 
thus k'g and X'^ are very different) will contribute strongly to 
the diffusion, while neurons for which kg ~ X\ will have little 
influence. 

If the population of neurons responds overall asymmet- 
rically to C, then one would still expect to observe lumpy 
response time histograms if the neurons are not too many. 
In the special case where X[ > AqVz (e.g., when C = 1 cor- 
responds to a sound and C = 0 corresponds to no sound, 
and all neurons are excited by sound), then all action poten- 
tials will send the diffusion upwards and the only way a NO 
decision may be taken is by the diffusion drifting toward 
the negative threshold. In this case the NO response time 
histogram will be sparse. However, the analysis of this case 
becomes considerably more tedious than the analysis presented 
in Section 5. 

If, on the other hand, the neurons are symmetrical w.r. to C, 
i.e., some neurons respond preferentially to C = 0 and some to 
C = 1, then the drift term is zero because drift terms for different 
neurons are equal and opposite and thus cancel each other. In this 
case the response time histograms are neither lumpy, nor sparse 
(see Figure 8). 

7. FIRING RATES AND RESPONSE TIMES 

It is intuitive that response times will be lower when the maxi- 
mum firing rate Xi is much larger than the minimum firing rate 
Xq. Consider Equation 8. Decisions based on zero action poten- 
tials take tq/S = xqKXq — Xi). Since tq is negative, these deci- 
sions are quicker when X\ >> Xq, which makes the magnitude 
of the denominator larger. Similarly, when fc > 0, decisions take 
are quicker when Xi >> Xq since the numerator increases log- 
arithmically, while the denominator increases linearly. Figure 9 
shows the behavior of RT vs Amax when Xmin = 1 for the sym- 
metric case of two neurons described in Section6, where for one 
neuron Xi = Xmax and Xq = Amin and vice-versa for the other 
neuron. 



8. METHODS 

The optimal Bayesian observer based on SPRT was implemented 
in Matlab using the Bernoulli method described in Section 3 
with Af = 1 ms and other parameters as shown in Table 1, unless 
otherwise specified in the text. 

9. DISCUSSION AND CONCLUSIONS 

The analysis I presented predicts that when binary decisions are 
computed by a mechanism involving a single neuron, one will 
observe sparse response time histograms for NO responses, and 
lumpy response time histograms for YES responses. 
A number of additional observations are possible: 

1. The histograms will still look highly structured when few 
neurons are involved, provided that the YES/NO decision is 
perceptually asymmetric (e.g., detecting the presence vs the 
absence of a sound, a vibration, or a light spot). 

2. The response time histogram wiU look like a log-normal dis- 
tribution when the neurons involved respond symmetrically 
to C. 

3. Even if the timing of the action potentials produced by the 
input neuron(s) are unpredictable and uncorrelated, as mod- 
eled by Poisson statistics, the timing of action potentials 
produced by the neurons that compute a decision is highly 
structured; it is tightly quantized in the case of NO decisions. 
Conceivably, this fact will enable additional computations 
where synchronization between action potentials is required. 

4. If quantized/structured response time histograms are observed 
in an experiment, it may be possible not only to predict that a 
small number of neurons is involved, but also to estimate their 
firing rate. 

5. The predictions I make are not affected by the real statistics of 
the neuron's response. They are a consequence of using SPRT 
and the assumption that the neuron's statistics is Poisson. It is 
impossible to decide whether a short sequence of action poten- 
tials is governed by Poisson or other statistics. I suspect that 
any Poisson-like assumption will lead to the same qualitative 
prediction. 

Observing sparse response time histograms experimentally is 
difficult for a number of reasons: 

1 . In most systems the number of neurons involved in the com- 
putation is more than a handful. 

2. If the task is symmetric (e.g., discriminate between a red and 
a green light), then there will be neurons tuned to both condi- 
tions making both YES and NO histograms continuous, rather 
than sparse. 

3. Motor response and neural propagation delays are them- 
selves a random variable which, if summed to the perceptual 
response time, will blur away the sparse nature of the NO 
histograms. 

I believe that one should be able to design single-neuron asym- 
metric preparations. For example, by stimulating a single gan- 
glion cell in the retina with a small light dot, or stimulating a 
single tactile receptor in the skin of the back, where receptors 
are sparse. If such experiments prove to be possible, response 
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time histograms would provide a wealth of information on the 
mechanisms involved in decision, including the firing rate of the 
neuron involved. 
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